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Abstract 

The two-body knock-out reaction 4 He(e, e'd)d is calculated at various momentum transfers. 
The full four-nucleon dynamics is taken into account microscopically both in the initial and the 
final states. As NN interaction the central MT-I/III potential is used. The calculation shows a 
strong reduction of the coincidence cross section due to the final state interaction. Nonetheless the 
theoretical results exhibit a considerable overestimation of the experimental cross section at lower 
momentum transfer. Comparisons with other, less complete, calculations suggest that consideration 
of a more realistic ground state might not be sufficient for a good agreement with experiment, rather 
a more realistic final state interaction could play an essential role. 
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I. INTRODUCTION 



Electron induced two-body knock-out reactions are considered to be an important tool 
to investigate nucleon-nucleon (NN) correlations in nuclei. In comparison to one-nucleon 
knock-out processes two-nucleon emission reactions give more detailed information on NN 
dynamics, while many important details are already integrated out in the former case. 
Hence quite a number of experimental and theoretical studies have been devoted to (e, e'NN) 
reactions in order to study NN correlations (see e.g. [1]). However, it is often not very easy 
to obtain a clear picture of the two-body correlations, in fact one would need theoretical 
calculations where all relevant effects contributing to the observable under investigation 
are taken into account. For this reason in many theoretical works also effects from final 
state interaction (FSI), as well as meson exchange and A currents have been considered. 
Microscopic calculations of two-body knock-out reactions with the proper final continuum 
state have only been carried out for two- and three-nucleon systems (see e.g. In 
more complex nuclei such exact and consistent studies are still missing. 

The aim of the present paper is as follows. We want to consider the full FSI in a two-body 
emission reaction microscopically and at the same time go beyond the above mentioned two- 
and three-nucleon systems. To this end we consider the 4 He(e, e'd)d reaction. The choice of 
this particular reaction is based on three different considerations. The first is that, differently 
from the two- and three-body systems, 4 He has some characteristics (binding energy per 
nucleon, density) rather similar to those of heavier nuclei. Secondly, the (e, e'd) reaction 
has been suggested as a particularly useful tool for investigating short range correlations 
(for a brief summary see |^|). In fact the 4 He(e, e'd)d reaction was among the first NIKHEF 
experiments dedicated to the study of NN correlations 5(. The third reason is that in ref. 
[f| large effects of FSI were found in a two-deuteron cluster model. 

In the present paper we perform a rigorous calculation of the FSI of the four-nucleon 
system, where the nucleons interact via an NN potential. We employ an integral transform 
method as outlined in ref. [7] . Particularly suited for such kind of calculations is the Lorentz 
Integral Transform (LIT), which was proposed in ref. [8(. In fact the LIT has alr eady 
been used for the calculation of various exclusive reactions: d(e,e'p)n jjj, 4 He(7,p) 3 H jlQ| . 
4 He( 7 ,n) 3 He Q and 4 He(e, e'pfH. Q. 

Our paper is organized as follows. In sect. 2 we describe the formalism of the 4 He(e, e'd)d 
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cross section. A discussion of results and a short conclusion are given in sect. 3. 



II. FORMALISM 

The longitudinal part of the 4 He(e, e'd)d coincidence cross section is given by 



de k >dVL e ,dVL d " 1V1 2(1 - (q/2p d ) cos(0 d )) q 4 



d 5 a L M dPd q^ 

°MTT, 7 /o \ 77Ty\ a F L{^,q,O d ). (1) 



We do not consider other cross section parts, where the transverse current is involved. The 
reason is that we mtelld to compare our results to the data o f r e, g lu Q it was shown 
that for the kinematics of that experiment the transverse current contributions to the cross 
section are small. In eq. (JTJ ey and VL e i denote energy and solid angle of the scattered 
electron, and om is the Mott cross section. Energy and momentum transfer to the nuclear 
system are denoted by u and q = gq, q 2 is the squared four-momentum transfer, 9 d denotes 
the angle between p d = p d p d , the outgoing deuteron momentum, and q, while M d is the 
deuteron mass. Note that for given values of u, q, and p d there is a unique deuteron knock- 
out angle 9 d . The quantity Fi,(uj,q,9d) is the longitidinal structure function, defined as 

F L (q,u>,O d ) = £ |<<W^) lp(q)| * a }f ■ (2) 

M,M' 

Here ^m M , is the internal continuum final state of the minus type pertaining to the d-d 
channel with a given relative momentum of the final d-d pair (denoted below by k) and 
deuteron spin projections M and M' . The quantity E djd denotes the excitation energy of 
the four-nucleon system, 

k 2 

E d ,d = £d,d - 2E d + E a with e d:d = — , (3) 

M d 

where E d is the deuteron binding energy and E a is the 4 He binding energy. The four-body 
ground state is denoted by ty a , and G P E is the proton electric form factor. As nuclear charge 
operator p we take 

4 1 + r 3 

Here denotes the third component of the j-th nucleon isospin and rj represents the 
position of the j-th nucleon with respect to the center of mass of the four-body system. 
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In our calculation we do not make use of the continuum wave function $?m M i, but instead 
determine FL(q,u>,8 d ) by means of the integral transform method for exclusive reactions 
j]| with the Lorentz kernel j^, Q]. In this approach one avoids to treat a continuum state 
problem, one works instead with a much easier bound-state-like problem. Our calcu 
carried out in complete analogy to the 4 He(7,p) 3 H and 4 He(e, e'p) 3 H calculations of 
thus in the following we give only a very brief summary of the method. 

The starting point of the calculation are the transition matrix elements 



at ion is 
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111, 



TMM'(Ed t d) 

They can be divided into a Born term, 
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rpBorn ( 771 
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d,d) 



MM 



and a FSI dependent term, 



rpFSI ( pi 
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(5) 
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(7) 



E d ,d + ie - H 

In the above two equations <Pmm' ^ s a product of the internal wave functions of the two 

deuterons, with spin projections M and M', and of the Coulomb function for their relative 

motion with a given large-distance momentum. The relative motion occurs in the average 

d-d Coulomb potential. The potential V is the sum of all interactions between nucleons 

belonging to different deuterons with the average d-d Coulomb potential being subtracted, 

H denotes the full Hamiltonian of the four-nucleon system, and A is an antisymmetrization 

operator explicitly given in ref. 

The calculation of the Born term is rather unproblematic, while the FSI term is calculated 

with the integral transform method. To this end first the following identity is used 

f°° F MM i(E) 
! E - E d>d + ie 



rpFSI 

1 MM 



d,d) 



E 



-dE 
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where E t h is the break-up threshold energy of 4 He in the isospin T = channel and 



MM 



dvU- MM ,{E d4 ) VA y„(E u )) (y u {E v )\p\V a )5(E - E u ) 



(9) 



The function Fmm 1 contains information on all the eigenstates ^ v for the whole eigenvalue 
spectrum of H. In the LIT method it is obtained by its Lorentz integral transform 

Fmm 1 (E) 



L [F M m>] {a) 



E - (E - a R f + aj 



dE 



* 2 (<t) I #l(<7 



(10) 



where 

= (H — <j)- l p |* a ) , tt 2 (cr) = (H- a)- l AV\r M KAE d , d )) (11) 

and a = cr R + iaj. Equation (fTUj) shows that L [F MM r] (a) can be calculated without explicit 
knowledge of F MM i, provided that one solves the two equations 



(H-a) ^) = p\* a ) , (12) 
(H-a) * 2 ) = AV\r MM >(Ed,d)) ■ (13) 



The quantities and \&2 have finite norms and thus only bound state techniques are 
required to obtain the solutions of eqs. (JT2J) and (JT3J). 

We use expansions over a basis set of localized functions consisting of correlated hyper- 
spherical harmonics (CHH) multiplied by hyperradial functions, which lead to rather large 
Hamiltonian matrices. Instead of using a time consuming inversion method we directly 
evaluate the scalar products in (|1(J|) with the Lanczos technique as explained in ref. 

After having calculated L[Fmm'](c) one obtains the function Fmm'(E), and t 
TMM'{Ed,d), via the inversion of the LIT. We perform the inversion as described in 
(for other inversion methods see ^J). 
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III. RESULTS 

In our calculation we use the semi-realistic MT-I/III potential [16] as NN interaction. 
Below pion threshold it leads to rather good descriptions of the NN s-wave phase shifts 
3 S\ and 1 Sq. The Coulomb interaction is considered in addition. As already mentioned, the 
ground state of 4 He as well as ^ of eqs. (fE^ and (pH?j) are calculated using the CHH expansion 
method. In order to speed up the convergence, state independent correlations are introduced 
as in 0] (our 4 He wave function is identical with the CHH ^ a of Q]). The deuteron 
ground state wave function is determined by a numerical solution of the radial Schrodinger 
equation. As proton electric form factor we take the usual dipole parametrization. The 
transition matrix elements (JBI) and (|IJ) are calculated in the form of partial wave expansions. 
For the Born term we take into account multipoles up to order L = 20, while for the FSI 
term we include multipoles up to L = 5. We checked that with such an expansion a sufficent 
convergence is reached. 
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In this work we consider the 4 He(e, e'd) reactions measured in the above mentioned 

n 

NIKHEF experiment |5J. This corresponds to the following kinematical settings: relative 
energy of the two final deuterons e^d — 35 MeV and four-momentum transfers q 2 = 1.75, 
2.49, 3.36 and 4.79 fm" 2 . 

We start our discussion considering the effect of the Coulomb FSI. In fig. 1 results are 
shown for the highest of the four considered momentum transfers. As one may expect the 
effect is very small and only visible in the cross section minimum around 9d = 50 degrees. 
Similarly small effects of the Coulomb FSI are found for the other three momentum transfers. 
Also shown in fig. 1 is the PWIA result of [4;J, where a harmonic oscillator (HO) s-state wave 
function is taken for the 4 He ground state. Although such a ground state model is rather 
different from ours, one finds rather similar results in both calculations for the cross section 
minimum, while at forward and backward angles the HO ground state leads to an increase 
of the cross section by about a factor of two. 

In fig. 2 we illustrate the effect of the full FSI for all four considered momentum transfers. 
Though the positions of the minima at about 50 degrees are hardly changed, it is readily 
evident that FSI is very important for all four considered momentum transfers. For the lower 
two q-values one obtains reductions of strength by about a factor of 5, while the shapes 
of the angular distribution remain almost unchanged. On the contrary for the two higher 
momentum transfers also the shape is affected. One finds stronger decreases of Fl only in 
the forward region and for deuteron angles beyond 100 degrees and accordingly much less 
pronounced minima. 

In fig. 3 we show the differential cross section resulting from the calculated Fl at q 2 , = 4.79 
fm -2 . Different from the PWIA case the calculation with inclusion of FSI does not exhibit a 
minimum, but shows a rather constant fall-off with increasing deuteron knock-out angle. As 
already seen in fig. 2 the cross section is significantly reduced by FSI for forward directions 
and beyond 100 degrees. In comparison to the experimental cross section at $d — 15 degrees 
one finds for the PWIA an overestimation by about a factor of 2.6 and an underestimation 
by about a factor of 1.7 for the full calculation. 

In fig. 4 we make a comparison of our results with experimental data also for the other 
considered momentum transfers. One sees that the rather strong reductions of the cross 
sections due to FSI is by far not sufficient for an agreement with experiment at the lower 
momentum transfers. In fact the disagreement with experiment becomes more and more 



6 



pronounced with decreasing momentum transfer. For the lowest q-value the experimental 
cross section is overestimated by almost an order of magnitude. We mention once again that 
an inclusion of the here not considered transverse current contributions should not lead to 
significant changes. In fact in ref. [(| it was shown that the one-body current increases the 
PWIA result by less than 20 % for the lower three momentum transfers and by about 25 % 
for q 2 = 4.79 fm~ 2 . Additional two-body currents are expected to be small because of the 
isoscalar nature of the 4 He(e, e'd)d reaction. In fig. 4 we also show theoretical results from 

n 

|6J. A considerably smaller discrepancy between theory and experiment is found than in our 
calculation. It should be pointed out that in jf| FSI is not calculated microscopically as in 
our case, but modelled via a central phenomenological potential between the two outgoing 



deuterons. As already mentioned the 4 He ground state of |[6j is a simple harmonic oscillator 
wave function. 

In order to better understand the origin of the differences between the two calculations 

we show in fig. 5 the corresponding PWIA results. One notices that a large part of the 

differences of fig. 4 are already found for the PWIA results and thus are caused by the 

different ground state wavefunctions. In the figure we also illustrate the PWIA calculation 

of Morita taken from 161, where correlations are introduced in the 4 He ground state via the 
□ 

ATMS method DJ. The ATMS and HO results are rather similar and considerably lower 
than ours below q 2 = 2.5 fm -2 . The different results show that it is not very likely that a 
fully realistic calculation of the He ground state could close the gap to experiment at lower 
momentum transfer. Most probably one needs a more realistic description of the FSI for 
the four-nucleon final state. It should be noted that the FSI effects of the present work 
and those of Q| are not very different, in fact in both cases one finds similar reductions of 
the corresponding PWIA cross sections. However, we want to emphasize once again that 
there are principal differences for the treatment of FSI in both calculations. On the other 
hand both FSI calculations have in common that only central potentials between nucleons 
or respectively deuterons are taken into account. Probably a consideration of tensor and 
other realistic force terms could change the picture significantly. 

We summarize our work as follows. We have calculated the 4 He(e, e'd)d reaction tak- 
ing into account the full four-nucleon dynamics in inital and final states. Such two-body 
knock-out reactions are considered as an excellent tool to determine two-body ground state 
correlations. In fact the comparison of our PWIA results with corresponding results from 
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other calculations show a non negligible sensitivity of the cross section to ground state cor- 
relations. On the other hand we find that FSI effects are even more important leading 
to a strong reduction of the cross section. We thus confirm similar results from ref. jj 
where, however, FSI was considered in a d-d cluster model only and not as in our case mi- 
croscopically via a NN interaction. Compared to experimental data our theoretical results 
show deviations up to about 50 % in the range 3.3 fm -2 < g 2 < 4.8 fm~ 2 , but exhibit a 
considerable overestimation of the experimental cross section at lower momentum transfers. 
As already pointed out a more realistic 4 He ground state will most likely not be sufficient 
for a satisfying improvement, whereas a more realistic nuclear interaction also for the final 
four-nucleon state might lead to an agreement of theoretical and experimental results. 
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FIG. 1: Angular distribution of F L at e AA = 35 MeV and <? 2 = 4.79 fm -2 without any FSI (dashed 
curve) and with inclusion of Coulomb-FSI (full curve); also shown the PWIA result of 4] (dotted 
curve) . 
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FIG. 2: Angular distribution of Fl with (full curve) and without (dashed curve) FSI contribution 
at ed,d = 35 MeV and momentum transfers as indicated in the figure. 



11 



PWIA 
PWIA+FSI 
NIKHEF 



20 40 60 80 100 120 140 160 180 

d (deg) 



FIG. 3: Angular distribution of differential cross section at tdd = 35 MeV and qi = 4.79 fm 



with (full curve) and without (dashed curve) FSI; also shown experimental result (dot) from 



0. 



12 



10 



-5 



% lO" 6 



10 



-7 



10" 



10" 



1.5 



NIKHEF • 



PWIA 
PWIA+FSI 
Cluster Model 



2.5 



3.5 

q 2 u [fm" 2 ] 



4.5 



FIG. 4: Differential cross section at e^d = 35 MeV and missing deuteron momentum p m = q— Pd of 
|p m | = 125 MeV/c as function^of with (full curve) and without FSI (dashed curve) in comparison 



to experimental data (dots) [5j; also shown result from |(| with an HO 4 He ground state and FSI 
in a d-d cluster model (dotted curve). 



13 




14 



